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Among  the  items  considered  are  target  characteristics  that  affect 
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coherent  and  incoherent  illumination  with  respect  to  resolution. 

In  addition,  this  paper  considers  the  possibility  of  improving  image 
quality,  particularly  resolution,  with  the  aid  of  superresolution  tech¬ 
niques  and  apodization.  It  also  estimates  the  possible  reduction  in  the 
number  of  elements  in  a  large  antenna  array  by  the  use  of  an  optimum 
array  thinning  distribution  of  the  elements. 
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EXECUTIVE  SUMMARY 


This  paper  explores  the  theoretical  limits  on  performance 
of  radar  Imaging*  systems  whose  purpose  is  to  recognize  by 
shape  and,  as  far  as  possible,  to  identify  military  targets  on 
the  ground.  The  analysis  stresses  the  3.2  mm  wavelength  case, 
and  tanks  are  assumed  to  be  typical  targets  of  interest. 

The  major  conclusion  is  that  millimeter-wave  radars  cannot 
provide  imagery*  of  sufficient  quality  to  identify  or  even 
recognize  tank-sized  targets  at  tactically  interesting  ranges 
without  making  use  of  antennas  that  are  too  large  to  be  con¬ 
sidered  feasible.  A  possible  exception  is  the  synthetic  aper¬ 
ture  radar,  which,  however,  requires  an  airborne  platform,  and 
for  which  the  image  quality  nill  be  vulnerable  and  highly  sen¬ 
sitive  to  unpredlcted  target  motion. 

In  addition  to  imaging  system  parameters,  the  investigation 
considered  target  characteristics  that  affect  image  quality. 

It  also  considered  what  Image  quality  is  needed  for  target 
recognition  and  identification. 

When  the  surface  of  an  extended  target  is  entirely  smooth, 
imaging  is  not  possible.  Ideally,  for  imaging  the  surface 


In  the  context  of  this  paper  "image"  is  to  be  understood  in 
the  narrow  sense  of  ordinary  usage.  That  is,  the  term  refers 
here  to  the  kind  of  geometrical  representation  that  the  human 
eye  ordinarily  views  and  interprets.  Intentionally  excluded 
are  more  general  means  of  identification  that  are  often  pro¬ 
posed,  such  as,  for  example,  the  correlation  of  various  radar 
measurements  with  known  target  signatures. 
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should  be  uniformly  rough  relative  to  the  wavelength.  That  is, 
it  should  be  covered  densely  with  irregularities,  on  each  of 
which  there  is  a  distinct  specular  point  of  reflection.  Al¬ 
though  the  irregularities  and  the  spacing  between  them  must  be 
smaller  than  a  resolution  element,  at  each  specular  point  the 
principle  radii  of  curvature  of  the  target  surface  should, 
nevertheless,  be  larger  than  about  1.5  wavelengths. 

The  number  of  distinct  scattering  centers  that  can  exist 
on  a  target  is  limited  by  the  wavelength.  When  the  resolution 
of  an  imaging  system  is  sufficiently  fine,  the  image  provided 
by  the  return  from  those  scattering  centers  will  not  appear 
filled  but,  rather,  will  appear  to  consist  of  separated  points. 
This  suggests  that  once  such  a  resolution  has  been  reached 
additional  fineness  may  be  of  little  benefit. 

Thus,  the  value  of  achieving  fine  resolution  will  depend 
upon  the  target.  For  example,  aircraft  tend  to  be  relatively 
smooth  and  therefore  may  be  difficult  to  recognize  even  with 
the  fine  resolution  for  which  millimeter-wave  systems  have  a 
potential  capability.  On  the  other  hand,  existing  scattering 
data  indicate  that  ground  vehicles  such  as  tanks  will  provide 
improved  imagery  with  the  resolution  improvement  that  milli¬ 
meter-wave  radar  can  achieve. 

As  long  as  the  signal-to-noise  ratio  is  good,  the  resolu¬ 
tion  needed  for  identifying  tank-sized  targets  is  such  as  to 
provide  about  6  line  pairs  across  the  smallest  target  dimension. 
For  the  case  of  primary  interest  here,  this  is  equivalent  to  a 
requirement  that  the  imaging  system  be  able  to  detect  modula¬ 
tion  at  spatial  frequencies  at  least  as  high  as  3  cycles/m, 
which  Is  equivalent  to  a  minimum  spatial  period  at  least  as 
small  as  0.33  m.  The  requirement  for  recognition  is  not  sig¬ 
nificantly  lighter,  since  experiments  show  that  for  that  pur¬ 
pose  about  4  line  pairs  across  the  target's  smallest  dimension 
are  needed. 


A  review  of  what  is  known  in  general  about  imaging  indi¬ 
cates  that  predictions  sometimes  made  concerning  image  quality, 
based  on  classical  definitions  of  the  optical  resolution  limit, 
may  be  too  optimistic.  Two  reasons  for  this  are  evident:  fail¬ 
ure  to  take  into  account  the  effect  of  the  signal-to-noise  ratio 
and  failure  to  recognize  the  difference  in  effect  between  co¬ 
herent  and  incoherent  illumination. 

In  this  connection  it  is  important  to  distinguish  between 
extended  targets  and  point  targets.  The  classical  resolution 
limits  refer  to  point  targets.  They  are  valid  for  extended 
targets  only  when  the  illumination  is  incoherent  and  the  signal- 
to-noise  ratio  is  large,  i.e.,  of  the  order  of  20  dB  or  more. 

In  characterizing  the  image  quality  for  an  extended  target,  the 
modulation  transfer  function  is  more  pertinent  than  resolution. 

This  paper  also  considers  optimization  methods  that  might 
extend  the  limits  on  image  quality  or  contribute  to  the  feasi¬ 
bility  of  an  imaging  system.  Some  attention  was  given  to  the 
possibility  of  improving  resolving  power  by  means  of  super- 
resolution  techniques  and  to  the  question  of  how  many  elements 
would  be  required  for  an  antenna  array  if  nonuniform  spacing 
of  the  elements  were  permitted. 

Investigation  disclosed  that  superresolution  techniques 
will  not  improve  the  image  quality  of  an  extended  target,  al¬ 
though  apodization  (aperture  weighting)  may  reduce  the  effect 
of  diffraction  noise  by  distributing  it  more  evenly.  This  is 
similar  in  effect  to  the  standard  technique  of  eliminating  the 
speckle  phenomenon  associated  with  coherent  illumination  by 
dithering  the  target  or  performing  an  equivalent  operation  on 
some  part  of  the  Imaging  system. 

A  calculation  demonstrates,  on  the  other  hand,  that  it  is 
feasible  w.  .ad*  e  the  number  of  elements  required  in  an  array 
of  sensors  for  a  staring  system  (i.e.,  one  that  does  not  scan 


the  field  of  view)  by  distributing  the  sensors  nonuniformly 
over  the  aperture.  If  the  field  of  view  is  narrow,  e.g., 
restricted  to  the  angle  subtended  by  a  single  target,  the 
reduction  from  the  number  of  sensors  that  would  be  required  if 
the  distribution  were  uniform  can  be  greater  than  50  percent. 


I.  INTRODUCTION 

A.  PRIMARY  OBJECTIVE  AND  RESULTS 

This  paper  is  a  review  of  the  limitations  on  performance 
that  are  inherent  in  radar  imaging  systems,  particularly  if 
their  purpose  is  to  identify  military  targets  on  the  ground. 

The  emphasis  is  on  3.2  mm  waves  used  against  tank-sized  targets. 

The  primary  conclusion  of  the  review  is  that  the  only  type 
of  millimeter-wave  system  with  a  reasonable  antenna  size  that 
might  accomplish  the  stated  purpose  is  the  synthetic  aperture 
radar.  With  more  conventional  antennas  the  required  imaging 
performance  is  not  feasible  at  tactically  useful  ranges. 

Even  at  a  range  of  only  1  km,  a  3.2  mm  radar  would  need  an 
effective  antenna  diameter  of  nearly  20  m  to  identify  a  tank. 

If  such  an  antenna  were  implemented  as  an  array  of  sensors  and 
the  illumination  were  restricted  to  a  field  of  view  covering 
not  much  more  than  a  single  target,  the  array  would  still  con¬ 
tain  a  minimum  of  nearly  17,000  elements. 

B.  SCOPE  AND  ANALYTICAL  APPROACH 

This  paper  will  examine  a  number  of  questions  keyed  to  the 
problem  of  identifying  a  tank  from  imagery  produced  by  a  3.2  mm 
radar  at  a  range  of  1  km.  The  reason  for  working  with  this  very 
short  range  is  that  results  based  on  It  may  be  easily  scaled  to 
values  corresponding  to  more  realistic  situations.  In  addition, 
the  results  will  act  as  a  conservative  estimate  or  bound  on  the 
feasibility  of  using  millimeter-wave  radar  imaging  systems  for 
the  intended  applications. 
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It  is  assumed  throughout  that  propagation  conditions  are 
ideal,  i.e.,  that  the  relative  index  of  refraction  of  the  at¬ 
mosphere  is  unity  and  remains  constant  over  space  and  time. 
Therefore,  losses  in  signal  strength  due  to  absorption  and  Mie 
scattering  still  remain  to  be  taken  into  account.  This  is 
also  true  of  phase  errors  due  to  variations  in  the  refractive 
index  and  phase  fluctuations  due  to  atmospheric  turbulence. 

At  the  relatively  short  ranges  envisioned  here,  spatial 
variations  of  the  index  of  refraction  should  not  be  significant. 
A  small  amount  of  turbulence  might  even  be  helpful  to  a  system 
that  relies  upon  coherent  illumination  because  it  would  intro¬ 
duce  phase  averaging,  which  would  tend  to  reduce  the  speckle 
effect  normally  associated  with  coherent  imagery. 

The  term  "image"  used  in  connection  with  radar  sometimes 
includes  within  Its  scope  nothing  more  than  a  set  of  signatures 
by  means  of  which  a  target  can  be  classified.  In  this  paper 
the  usage  Is  much  narrower.  Implying  that  the  final  output  of 
a  system  defined  as  imaging  is,  In  fact,  the  kind  of  geometri¬ 
cal  representation  that  the  human  eye  ordinarily  views  and 
interprets. 

Chapter  II  considers  those  characteristics  of  a  target 
that  make  it  a  suitable  candidate  for  surveillance  by  a  radar 
Imaging  system,  as  well  as  the  image  quality  needed  for  target 
Identification  or  recognition.  Chapter  III  discusses  the  sys¬ 
tem  parameters  that  affect  image  quality  and  estimates  the 
limiting  parameter  values  that  are  required  for  target  identi¬ 
fication  performance.  Chapter  IV  contains  a  brief  summary  of 
conclusions  resulting  from  the  observations  and  discussions  of 
Chapters  II  and  III. 

Appendix  A  is  a  tutorial  review  of  general  image-process¬ 
ing  theory.  It  is  included  as  a  convenient  reference  for  the 
basic  relations  and  estimates  presented  in  Chapter  III  and  as 
background  material  for  the  array-thinning  analysis  presented 
in  Appendix  B. 
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An  attempt  has  been  made  to  consider  all  avenues  for  re¬ 
ducing  system  requirements  and  for  improving  the  fundamental 
parameters  used  to  measure  image  quality.  Those  avenues  include, 
in  the  first  instance,  minimizing  the  number  of  elements  needed 
in  an  array  of  sensors  and,  in  the  second,  refining  resolution 
by  means  of  superresolution  techniques.  No  attempt  has  been 
made,  however,  to  analyze  or  discuss  processing  techniques 
that  aim  at  enhancing  the  effectiveness  of  the  image  display, 
i.e.,  techniques,  such  as  selective  area  averaging  or  direc¬ 
tional  filtering,  that  may  properly  be  classified  under  the 
headings  of  optical  data  reduction  or  image  enhancement. 


II.  TARGET  CHARACTERISTICS 

A.  THE  NATURE  OF  RADAR  SCATTERING  AND  IMAGE  FORMATION 

If  the  electromagnetic  radiation  due  to  a  monochromatic 
radar  source  is  reflected  from  a  large,  smooth,  plane  mirror, 
the  reflected  field  will  appear,  in  all  respects,  to  come  from 
the  geometrical  image  of  the  source  relative  to  the  mirror. 
Therefore,  the  radar  return  will  contain  no  information  about 
the  shape  of  the  mirror  except,  perhaps,  what  can  be  inferred 
indirectly  from  the  fact  that  the  image  of  the  source  is  un¬ 
distorted. 

If  the  mirror  is  curved — for  example,  if  it  is  the  surface 
of  a  large  sphere — the  return  will  also  appear  to  come  from  an 
image  of  the  source.  However,  in  this  case  the  image  will  be 
distorted  and  blurred.  The  effect  of  the  change  in  shape  of 
the  mirror  surface  is  to  introduce  optical  aberrations.  Never¬ 
theless,  It  Is  still  the  source  that  Is  imaged,  not  the  reflec¬ 
tor  surface. 

On  the  other  hand,  if  the  reflector  Is  not  smooth  but  is 
covered  uniformly  with  small  irregularities,  energy  will  appear 
to  be  scattered  separately  from  points  evenly  distributed  over 
the  surface  wherever  It  is  illuminated.  Macroscopically ,  the 
return  will  appear  to  come  from  an  extended  region  that  is  geo¬ 
metrically  similar  to  the  illuminated  part  of  the  reflector. 

When  the  radar  wavelength  is  small,  geometrical  optics  pro¬ 
vides  a  valid  theoretical  model  of  electromagnetic  scattering. 
This  Is  true  even  when  the  reflector  surface  is  rough,  as  long 
as  the  principal  radii  of  curvature  of  most  points  on  the  surface 
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are  larger  than  about  1.5  wavelengths.  At  the  other  extreme, 
microscopic  irregularities,  i.e.,  those  smaller  than  an  eighth 
of  a  wavelength,  have  a  negligible  effect  on  scattering. 

Thus,  it  is  reasonable  to  expect  that  roughness  may  be 
characterized  by  the  existence  of  a  large  number  of  specular 
points,  relative  to  an  Imaging  system,  due  to  the  irregularities 
distributed  over  the  reflector  surface.  In  fact,  statistical 
models  of  quasi-dif fuse  scattering,  in  good  agreement  with  ex-  . 
periment,  have  been  obtained  by  assuming  that  a  rough  surface 
consists  of  randomly  oriented  surface  elements,  each  of  which 
in  accordance  with  its  own  orientation  reflects  in  a  specular 
direction  relative  to  the  optical  rays  incident  on  it  (cf.  Ref.  1). 

At  each  of  the  many  specular  points  on  a  rough  reflector 
surface,  a  pencil  of  rays  from  the  incident  field  generates  a 
pencil  of  reflected  rays  that  pass  through  the  entrance  pupil 
of  the  imaging  system.  A  perfect  system  will  process  the  re¬ 
sulting  phase  distribution  over  its  aperture  to  form  an  image 
of  the  caustic*  associated  with  each  reflected  pencil  of  rays. 

If  rays  In  the  pencil  are  those  associated  with  a  spheri¬ 
cal  wave,  the  caustic  Is  a  point.  In  that  case,  the  imaging 
system  will  produce  the  diffraction-limited  image  of  a  point. 

Tf  the  caustic  Is  not  totally  degenerate  but  consists  of  lines 
or  surfaces,  the  Image  may  still  appear  to  correspond  to  a  point 
(of  best  focus);  however,  it  will  also  appear  to  have  aberra¬ 
tions,  the  nature  of  which  will  depend  upon  the  caustic  geometry. 

f - 

A  caustic  is  the  envelope  formed  by  the  two-parameter  family 
of  rays  that  are  orthogonal  to  a  wavefront  surface.  The  caus¬ 
tic  is  mathematically  equivalent  to  the  locus  of  the  principle 
centers  of  curvature  of  the  wavefront  surface  and  Is,  itself, 
usually  a  surface  consisting  of  two  branches  Joined  along  a 
singular  curve  called  the  edge  of  regression.  When  the  wave- 
front  has  appropriate  symmetry  the  caustic  may  degenerate  into 
a  surface  and  a  line,  two  lines,  or  even  a  single  point,  which 
would  then  be  the  (perfect)  focal  point  of  the  rays. 
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When  the  F-stop  of  the  imaging  system  is  large,  aberrations 
do  not  greatly  affect  resolution,  at  least  for  target  resolution 
elements  near  the  center  of  the  field  of  view.  Thus,  when  the 
return  is  backscattered ,  the  fact  that  images  of  resolution  ele¬ 
ments  in  the  target  surface  are  actually  images  of  caustics 
should  not,  in  itself,  cause  much  deterioration  of  the  resolv¬ 
ing  power  of  the  system.  If  the  illumination  is  bistatic, 
however,  some  deterioration  may  occur  because  of  coma  and  greatly 
increased  astigmatism. 

For  purposes  of  good  imagery  a  target  surface  should  be 
uniformly  rough;  that  is,  at  any  aspect  the  distribution  of 
specular  points  should  be  so  dense  that  there  will  be  no  regions 
more  than  a  resolution  element  or  so  wide  that  are  devoid  of 
specular  points.  There  are  two  reasons  for  this  requirement. 

One  is  the  desirability  of  having  enough  point  images  to  fill  in 
the  geometrical  shapes  that  characterize  the  object.  The  other 
is  to  avoid  the  highly  aspect-dependent,  excessively  large  return 
that  may  be  produced  by  a  smooth  surface  patch,  causing  glint 
or  glare  which  could  disrupt  the  scene  or  even  strain  the 
system’s  dynamic  range  capacity. 

Other  properties  of  the  target  may  affect  image  quality, 
such  as  lossy  surface  materials  that  reduce  reflectivity. 
Reflectivity  may  also  be  highly  dependent  on  polarization, 
e.g.,  when  multiple  reflections  take  place  in  the  scattering 
process.  The  dependence  of  reflectivity  on  polarization  is 
most  pronounced  for  bistatic  illumination,  and  depolarization 
can  occur  in  that  case  even  in  the  absence  of  multiple  re¬ 
flection. 

The  intensity  of  a  wave  reflected  from  a  specular  point 
depends  not  only  on  the  reflectivity  due  to  material  properties 
of  the  reflector,  but  also  on  the  local  curvature  of  the  surface. 
This  geometric  dependence  is  a  consequence  of  the  fact  that  the 
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intensity*  of  a  propagating  wave  varies  in  accordance  with  the 
classical  inverse  square  law  of  geometrical  optics. 

The  generalized  version  of  the  law  is  that  the  intensity 
at  any  point  on  a  wavefront  is  inversely  proportional  to  the 
product  of  the  wavefront's  two  principal  radii  of  curvature 
at  that  point.  Each  of  these  radii  of  curvature  is  the  dis¬ 
tance  from  the  wavefront  along  the  ray  through  the  point  to 
the  ray's  point  of  tangency  on  one  of  the  two  branches  of  the 
associated  caustic.  The  law  is  also  equivalent  to  the  conser¬ 
vation  of  energy  flux  in  every  infinitesimal  pencil  of  rays; 
i.e.,  it  is  assumed  that  no  energy  can  leak  from  one  infinite¬ 
simal  pencil  of  rays  to  another  and  that  energy  is  conserved. 

In  the  case  in  which  the  Incident  field  is  a  plane  wave 
propagating  in  a  direction  normal  to  the  reflector  surface, 
each  principal  radius  of  curvature  of  the  wavefront  on  the 
backscattered  ray  at  the  point  of  reflection  is  one- half  of  a 
principal  radius  of  curvature  of  the  reflector  surface.  This 
leads  to  the  well-known  formula  for  the  radar  cross  section 
a  due  to  specular  reflection  from  a  large  smooth  body 
(cf.  Ref.  2),  namely, 

a  =  trp1p2  ,  (1) 

where  p^  and  p2  are  the  principal  radii  of  curvature  of  the 
body's  surface  at  the  specular  point.  The  factor  it  in  Eq.  1 
results  from  the  definition  of  radar  cross  section,  according 
to  which  the  normalized  intensity  of  the  scattered  field  is 
multiplied  by  4ir. 


The  term  "intensity"  as  used  here  has  the  meaning  that  is 
common  in  radar  terminology,  where  it  is  synonymous  with  flux 
density  measured  in  watts  per  square  meter,  rather  than  the 
meaning  that  is  common  in  optics,  where  it  refers  to  the 
strength  of  a  source  and  is  measured  in  watts  per  steradian. 
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The  calculation  is  more  complicated  if  the  scattering  is 
bistatic  or  the  incident  wave  is  nonplanar.  However,  in  that 
case  it  is  still  true  that  the  intensity  of  the  scattered  field 
is  proportional  to  the  product  of  the  principal  radii  of  curva¬ 
ture  of  the  reflector  surface  at  the  specular  point. 

B.  REQUIREMENTS  FOR  TARGET  RECOGNITION  AND  IDENTIFICATION 

A  pragmatic  way  to  determine  what  image  quality — in  par¬ 
ticular,  what  resolution  and  contrast — is  needed  for  the  recog¬ 
nition  and  identification  of  a  given  type  of  target  is  to  ana¬ 
lyze  the  target's  distinctive  visual  characteristics  pertain¬ 
ing  to  shape.  In  general,  the  answer  will  not  only  vary  with 
the  target  but  also  with  its  aspect. 

In  the  case  of  a  tank  viewed  broadside  the  most  distinctive 
features  are  probably  the  typical  parallelogram  shape  of  its 
tread,  the  trapezoidal  shapes  above  the  tread,  and  the  helmet- 
like  turret.  It  is  tempting  to  include  the  long  cylindrical 
shape  of  its  one  or  more  guns.  However,  some  evidence  exists 
that  gun  barrels  at  most  aspects  may  be  absent  from  millimeter- 
wave  imagery  of  tanks,*  probably  because  gun  barrel  surfaces 
are  relatively  smooth  in  this  wavelength  regime. 

One  characteristic  that  may  be  useful  for  identification 
purposes  is  the  number  of  a  tank's  bogey  wheels.  Thus,  a  stan¬ 
dard  for  the  resolution  needed  to  identify  tanks  might  be  based 


See  Ref.  3,  which  contains  a  number  of  figures  showing  3.2- 
millimeter-wave  as  well  as  optical  imagery  of  the  same  tar¬ 
gets,  including  some  that  portray  a  tank.  Whenever  the  tar¬ 
get  is  a  tank,  the  guns,  although  present  in  the  optical, 
are  missing  from  the  millimeter-wave  imagery.  Reference  4** 
contains  a  data  base  of  synthetic  aperture  radar  imagery  on 
tanks  and  trucks.  Measurements  at  35  GHz  were  made  to 
characterize  reflection  centers. 

«« 

No  classified  material  from  this  classified  reference  has 
been  used  in  this  paper. 


on  what  would  be  required  to  count  its  bogey  wheels.  That  reso¬ 
lution  would  be  sufficient  for  estimating  lengths  and  widths  of 
other  parts  of  the  tank  to  within  a  tolerance  of  the  order  of 
a  bogey  wheel  radius. 

Another  approach  to  estimating  the  image  quality  needed 
to  identify  or  recognize  tanks  is  to  apply  known  results  of 
psychophysical  experiments,  such  as  the  number  of  resolution 
elements  that  must  be  perceived  in  a  target  image  for  recog¬ 
nition  or  identification.  Reference  5  discusses  in  detail 
such  data  and  conditions  under  which  it  may  be  used  reliably. 

Reference  5  lists  an  empirically  derived  set  of  rules, 
due  to  J.  Johnson  (Ref.  6),  that  prescribe  the  number  of  reso¬ 
lution  line  pairs  across  the  minimum  target  dimension  that  are 
required  for  a  hierarchy  of  discrimination  levels.  In  increas¬ 
ing  order,  these  levels  are  designated:  detection,  orientation, 
recognition,  and  identification.  A  cautionary  is  added,  em¬ 
phasizing  that  in  applying  the  rules  care  must  be  taken  not 
to  neglect  contrast  and  its  dependence,  along  with  resolution, 
on  the  signal-to-noise  ratio. 

Johnson's  rule  for  target  recognition  is  that  about  4  line 
pairs  of  resolution  must  exist  within  the  target's  smallest 
dimension.  For  identification  his  rule  increases  that  number 
to  about  6  line  pairs.  The  minimum  dimension  of  a  tank  is  gen¬ 
erally  its  height,  which  may  be  on  the  order  of  2  meters.  For 
tank  identification  this  translates  into  a  requirement  for  a 
target  resolution  corresponding  to  a  minimum  spatial  period* 
of  about  0.5  m  to  about  0.33  m,  which  is  comparable  to  the 
radius  of  a  bogey  wheel.  Thus,  Johnson's  rule  appears  to  be 
consistent  with  the  estimate  obtained  by  means  of  what  has 
been  characterized  here  as  the  pragmatic  approach. 

* - 

See  Section  III-B-2-e  for  a  discussion  of  the  relationship 

between  spatial  fequency  and  resolution. 
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ill. 


IMAGING  SYSTEMS 


A.  TYPES  OF  SYSTEMS 

Basically  -.here  are  two  types  of  radar  systems  that  image 
extended  targets.  One  is  the  line  scanning  system,  which 
actively  scans  the  target  surface  with  a  focused  beam,  inte¬ 
grating  the  separate  return  from  each  resolution  element  se¬ 
quentially  along  a  raster.  The  other  is  the  staring  system, 
active  or  passive,  which  receives  and  processes  simultaneously 
the  return  from  many  resolution  elements  on  the  target  surface, 
however  the  target  may  be  illuminated. 

For  a  line  scanning  system,  resolution  in  the  elevation 
and  azimuth  dimensions  depends  upon  the  beamwidth  of  the  illu¬ 
minator  and/or  receiver  in  the  focal  region  and,  therefore,  upon 
the  size  of  the  effective  aperture(s)  associated  with  the  device. 
For  a  staring  system,  the  resolution  depends  upon  the  size  of 
the  aperture  at  the  entrance  pupil  associated  with  the  receiv¬ 
ing  sensor. 

Either  type  of  radar  system  may  provide  imagery  over  another 
dimension  by  ranging.  For  this  purpose  various  techniques  are 
available:  time  gating,  pulse  compression,  and  a  variety  of 
heterodyning  schemes  that  can  be  used  with  CW.  In  principle, 
it  is  possible  to  estimate  the  range  of  a  resolution  element 
by  means  of  focusing;  however,  at  telescopic  distances  the 
depth  of  field  of  an  imaging  system  will  usually  be  too  large 
to  make  this  approach  practical. 

Either  type  of  system  may  use  an  array  of  sources  for 
active  illumination  or  an  array  of  sensors  for  reception.  At 
least  in  the  millimeter-wave  regime,  staring  systems  may  have 
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fewer  problems  with  arrays  than  scanning  systems.  For  passive 
reception  the  array  elements  can  remain  co-phased  relative  to 
a  single  local  oscillator,*  so  that  phase  shifters  would  not 
be  necessary. 

Scanning  systems  that  scan  in  two  dimensions  may  suffer 
from  ambiguities  which,  as  observed  in  Ref.  5,  can  result 
from  the  periodic  nature  of  the  raster.  Removing  such  am¬ 
biguities,  e.g.,  by  means  of  spatial  filtering,  could  degrade 
resolution  or  contrast. 

Synthetic  aperture  radar  (SAR)  is  an  imaging  system  having 
the  ability  to  produce  fine-resolution  imagery  with  a  small- 
aperture  antenna.**  Most  of  its  applications,  such  as  terrain 
mapping,  have  been  those  involving  broad-area,  two-dimensional 
targets . 

SAR  has  a  fundamental  limitation:  it  can  only  provide  fine- 
resolution  imagery  in  the  range  and  azimuth  dimensions.  For 
this  reason  and  because  of  the  fact  that  the  real  antenna  must 
travel  along  a  path  at  least  as  long  as  the  synthetic  aperture 
that  it  generates  during  the  integration  time,  it  is  likely 
that  the  platform  must  be  airborne  in  any  application. 

The  most  difficult  problem  for  a  SAR  is  generally  the 
phase  error  caused  by  deviations  of  the  platform  from  its  pre¬ 
scribed  path.  If  the  target  is  not  stationary,  the  phase  error 
due  to  target  notion  creates  an  equally  difficult  problem. 
Sensitivity  to  these  kinds  of  error  scales  with  wavelength; 
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See  Appendix  A  for  image  processing  details  that  Justify 
this  statement. 

# 

As  it  is  resolving  images  in  azimuth,  it  operates  basically  as 
a  staring  system.  However,  since  SAR  scans  in  range,  it  may 
alsc  be  regarded  as  a  linear  scanning  array. 
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therefore ,  a  millimeter-wave  system  should  be  an  order  of  mag¬ 
nitude  more  sensitive  to  errors  generated  by  target  and  platform, 
motion  than  an  X-band  system. 

In  mitigation,  it  should  be  pointed  out,  however,  that 
since  for  the  same  resolution  the  required  aperture  length  is 
shorter  for  millimeter  waves ,  errors  that  increase  or  accumu¬ 
late  in  time  will  have  less  time  to  do  so.  In  fact,  the  effect 
of  such  errors  can  generally  be  reduced  by  increasing  the  radar 
platform  velocity. 

SAR  systems  ordinarily  correct  for  some  of  the  error  cre¬ 
ated  by  platform  motion  deviations,  to  the  extent  that  it  is 
possible  to  sense  such  deviations.  In  principle,  with  the  aid 
of  range  rate  measurements  similar  corrections  can  be  made  for 
errors  due  to  the  gross  translational  motion  of  a  target.  How¬ 
ever,  the  system  appears  to  be  at  the  mercy  of  other  kinds  of 
target  motion  such  as  pitch  and  yaw,  and  must  rely  upon  shorten¬ 
ing  the  integration  time,  as  required,  to  reduce  their  effect. 

B.  RELATIONS  BETWEEN  IMAGE  CHARACTERISTICS  AND  SYSTEM  PARAMETERS 

1 .  Fundamental  Parameters  of  Imaging  Systems 

All  imaging  systems  are  subject  to  the  same  basic  limita¬ 
tions  and  tradeoffs  among  the  parameters  that  define  image 
quality  and  coverage.  The  possible  values  of  these  parameters 
depend  upon  a  small  number  of  fundamental  system  constants. 

They  are:  the  wavelength  X,  the  limiting  aperture  (assumed  here 
to  be  circular  with  a  diameter  D),  the  signal  bandwidth  1/t, 
and  the  signal-to-noise  ratio  S/N. 

2.  Image  Characteristics 

a.  Depth  of  Field.  If  a  diffraction-limited  imaging  sys¬ 
tem  is  focused  on  a  point  located  on  the  optical  axis  at  a  dis¬ 
tance  z  from  the  entrance  pupil  plane,  there  will  be  an  inter¬ 
val,  called  the  depth  of  field,  about  the  point  where  all  other 
points  are  also  in  focus  within  some  defined  resolution  criterion. 
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According  to  Appendix  A,  the  imaging  system  will  be  focused  on 
a  target  resolution  element  that  is  in  a  plane  a  distance  R  from 
the  entrance  pupil  plane  as  long  as  R  is  in  the  interval  de¬ 
termined  by 


z 


2Xz 


D  +  2Xz 


£  R  £  z  + 


2Xz‘ 


-  2Xz 


(2) 


2 

when  z  <  D  /2X,  where  z  is  the  distance  between  the  entrance 

2 

pupil  and  the  object  plane.*  When  z  s  D  /2X ,  the  upper  limit 
of  the  interval  is  ®  . 

For  the  aperture  diameter  D  sufficiently  large  compared 
with  the  wavelength  X,  the  interval  defined  by  Eq.  2  is  approx¬ 
imately  the  same  as  that  de_fined  by 

2Xz2  .  „  .  _  .  2Xz2 

z  _  s.  R  &  z  +  ,  (3a) 

D  D 


while  for  large  z  it  is  approximately  the  interval  defined  by 

n2 

Ijj-  s  R  <  »  .  (3b) 

Although,  in  principle,  the  fact  that  the  depth  of  field  is 
finite  could  provide  a  method  of  ranging,  it  is  clear  from 
Eq.  3  that  for  millimeter  waves  and  target  ranges  of  1  km  or 
more  the  range  resolution  would  be  too  poor  for  the  method  to 
have  any  practical  value.  Figure  1,  which  shows  the  variation 
of  the  depth-of-field  limits,  depicting  them  as  functions  of 
D  for  z  *  1  km  and  X  *  0.0032  m,  Indicates  this  graphically. 


I - 

The  subscript  on  the  variable  z  used  in  Appendix  A  has  been 
dropped  here. 
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DEPTH-OF  HELD  LIMIT 


The  derivation  in  Appendix  A  of  Eq.  2,  and  therefore 
Eq .  3,  depends  upon  the  Rayleigh  criterion  for  diffraction- 
limited  focusing.  Actually,  the  image  intensity  drops  19 % 
from  its  value  at  the  center,  i.e.,  the  position  of  best  focus, 
to  its  value  at  either  end  of  the  depth  of  field  interval. 

This,  indeed,  is  one  criterion  used  to  define  depth  of  field. 

b.  Range  Resolution.  Pine  range  resolution  can  be 
achieved  with  a  short  pulse.  If  the  pulse  width  is  t,  the 
return  can  be  divided  into  time  intervals  and  processed  for 
range  bins,  each  of  which  corresponds  to  a  slice  of  width  ^ 
out  of  the  target,  where  c  is  the  velocity  of  light.  Thus, 
a  1-nsec  pulse  will  provide  range  resolution  of  0.15  m. 

This  implies  a  signal  bandwidth  1/t  of  at  least  1  GHz. 
Longer  frequency-modulated  pulses  with  the  same  bandwidth  can 
provide  the  same  resolution  when  processed  by  means  of  pulse 
compression . 

Care  must  be  taken  to  avoid  range  ambiguities  over  the 
target.  If  the  maximum  target  dimension  In  range  is  AR,  the 
pulse  repetition  rate  fp  must  satisfy  the  condition 


fp  < 


(4) 


Thus,  for  a  target  extending  10  m  in  the  range  direction,  ambi¬ 
guities  may  occur  if  the  pulse  repetition  rate  exceeds  15  MHz. 
Of  course,  if  absolute  range  to  the  target  is  desired,  AR  must 
be  replaced  by  R  in  Eq.  4. 

c.  Object  Plane  Resolution  of  Point  Targets.  The  reso¬ 
lution  limit  rQ  of  two  point  targets  is  the  minimum  separation 
between  their  positions  such  that  they  can  still  be  recognized 
as  distinct  targets.  According  to  the  classical  Rayleigh  defi¬ 
nition,  if  the  illumination  is  Incoherent  and  the  object  plane 
is  a  distance  R  from  the  entrance  pupil  aperture  plane, 
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(5) 


16 


According  to  Ref.  7,*  when  the  illumination  is  coherent  this 
resolution  limit  is  degraded  by  a  factor  of  1.^2;  i.e.,  in  that 
case  the  resolution  limit  r^  is  given  by 


r 


» 
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1.73X 

D 


R 


(6) 


The  difference  between  the  resolution  limit  for  the  case  of 
coherent  illumination  and  that  for  the  case  of  incoherent  illu¬ 
mination  is  that  the  combined  pattern  for  coherent  illumination 
results  from  the  addition  of  amplitude  distributions,  whereas 
the  combined  pattern  for  incoherent  illumination  results  from 
the  addition  of  intensity  distributions. 

In  Ref.  8,  R.K.  Luneberg  defines  the  resolution  limit 
for  two  equally  bright  point  targets  as  the  separation  at  which 
their  image  diffraction  patterns  combine  to  form  a  total  dif  ft  avia¬ 
tion  pattern  which  has  a  single  point  of  maximum  intensity,  but 
which  for  larger  separations  would  have  two  points  of  maximum 
intensity.  The  resolution  limit  defined  in  this  way  is  equal 
to  the  distance  between  the  inflection  points  about  the  central 
maximum  of  the  pattern  for  the  image  of  a  single  point  target. 

The  resolution  limit  for  the  coherent  case  is  equal  to  the 
distance  between  the  inflection  points  about  the  central  maximum 
of  the  amplitude  pattern  for  the  Image  of  a  single  point  target. 
For  the  Incoherent  case  it  is  equal  to  the  distance  between  the 
inflection  points  about  the  central  maximum  of  the  intensity 
pattern  for  the  image. 

For  the  case  of  Incoherent  illumination,  Luneberg's  defi¬ 
nition  leads  to  the  relation 
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No  classified  material  from  this  classified  reference  has 
been  used  In  this  paper. 
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For  the  case  cf  coherent  illumination,  it  leads  to  the  relation 


_  1.60X  n 

ro  "  D  R 


(8) 


Luneberg's  definition  of  the  resolution  limit  implies  an 
error-free  measurement  of  intensity.  Thus,  the  existence  of  a 
finite  S/N  degrades  the  limit. 

Otherwise,  S/N  determines  the  maximum  number  v  of  distin¬ 
guishable  levels  of  intensity,  i.e.,  the  dynamic  range,  in 
the  image.  Standard  results  from  information  theory  provide 
estimates  of  v,  which  depend  upon  the  noise  statistics. 

For  the  case  of  white  noise,  for  example,  it  is  well  known 
that  the  information  content  o:  a  signal  amplitude  is  j  log2 
(1  +  S/N)  bits.  Since  the  sign  of  the  amplitude  is  lost  when 
only  the  intensity  is  known,  a  measurement  of  intensity  will 
contain  one  bit  less  information.  It  follows  that 


v  =  |  V  1  +  S/N  (9) 

because  the  information  in  bits  provided  by  a  measurement  of 
Intensity  is  log2  v. 

In  order  to  take  into  account  the  effect  of  S/N  on  reso¬ 
lution,  it  would  seem  natural,  therefore,  to  reformulate  the 
definition  as  follows .  The  resolution  limit  for  two  equally 
bright  point  targets  is  the  separation  for  which  the  combined 
diffraction  pattern  cf  their  images  has  two  points  of  eaual 
maximum  intensity  that  differ  from  the  iptensity  at  a  local 
minimum  between  them  by  an  increment  that  is  1/v  times  the 
intensity  at  the  maxima.  This  would  make  the  Rayleigh  definition 
for  the  case  of  incoherent  illumination  correspond  to  an  S/N  of 
about  20  dB,  since  the  minimum  Intensity  is  about  80$  of  the 
maximum  in  that  case. 
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The  resolution  limit  for  3AR  (focused)  is  usually  stated 
as  1/2,  where  D  is  the  diameter  of  the  real  antenna.  This 
estimate  is  actually  the  distance  between  half-power  points  in 
the  effective  beam  of  the  synthetic  antenna  at  the  range  of 
the  resolution  element  being  imaged.  However,  it  agrees  very 
closely  with  the  definition  given  by  Luneberg  for  extended 
targets  in  the  coherent  case,  i.e.,  Eq .  11  in  Section  III-B-2-e 
with  D  replaced  by  the  synthetic  aperture  length. 

d.  Superresolution.  A  number  of  authors  have  proposed 
schemes  for  improving  resolution  by  means  of  phase  or  amplitude 
weighting  of  the  return  signal  distribution  over  the  entrance 
pupil  or  some  other  system  aperture  (cf.  Ref.  9).  The  idea,  gen 
erally,  has  been  to  reduce  the  diameter  of  the  central  spot  in 
the  Airy  disc  diffraction  pattern  of  a  point  target  image,  sub¬ 
ject  to  some  optimization  condition  such  as  maintaining  the 
greatest  possible  peak  intensity  or  the  maximum  energy  in  the 
central  spot.  The  peak  intensity  is  always  degraded  by  such  a 
process  (Ref.  8),  and  energy  is  either  absorbed  or  added  to  the 
outer  rings  of  the  diffraction  pattern. 

This  type  of  processing  is  similar  in  some  respects  to  mono 
pulse  in  radar.  In  fact,  one  technique,  due  to  J.  E.  Wilkins 
(Refs.  9,  10),  is  the  optical  analogue  of  monopulse.  Wilkins 
suggests  applying  a  180-degree  phase  shift  to  the  aperture  dis¬ 
tribution  ever  the  central  part  of  the  pupil  within  a  circle  of 
specified  radius,  leaving  the  rest  of  the  pupil  undisturbed. 

The  circular  region  to  which  the  phase  change  is  applied  corre¬ 
sponds  to  the  monopulse  sum  pattern,  and  the  undisturbed  annular 
ring  in  the  pupil  corresponds  to  the  difference  pattern.  The 
result  is  a  smaller  diameter  for  the  central  spot  in  the  Airy 
disc,  with  the  peak  Intensity  optimized  for  its  greatest  pos¬ 
sible  value. 
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Reference  9  describes  a  nur.be r  of  superresolution  tech¬ 
niques  that  affect  the  point  target  image  diffraction  pattern 
in  different  ways.  Few  such  techniques  have  actually  been 
applied,  however,  although  there  are  some  obvious  candidates  for 
application. 

The  improvement  of  the  resolution  of  double  stars  has  been 
an  objective  of  image  processing  by  means  of  aperture  weighting 
in  astronomy;  however,  in  practice  this  seems  to  have  been  con¬ 
fined  to  apodization,  a  process  for  which  the  goal  is  to  reduce 
the  intensity  of  the  outer  rings  (sidelobes)  of  the  Airy  disc. 
Apodization  improves  image  quality  by  reducing  diffraction  noise, 
but,  technically,  the  effect  is  to  enhance  the  S/N  rather  than  to 
improve  resolution. 

Luneberg  (Ref.  8)  derived  several  aperture-weighting 
functions  for  superresclution ,  each  optimizing  the  diffraction 
pattern  with  respect  to  a  different  criterion.  However,  he  also 
demonstrated  mathematically  that,  using  aperture  weighting,  it 
is  impossible  to  increase  the  image  spatial  frequency  bandwidth, 
which  is  limited  by  the  system  modulation  transfer  function 
(MTF).  The  implication  is  that  while  superresolution  is  possible, 
in  principle,  for  point  targets,  it  is  not  possible  for  extended 
targets . 

e.  Extended  Target  Resolution.  In  Ref.  8  Luneberg  also 
proposed  an  alternative  definition  of  the  resolution  limit,  suit¬ 
able  for  an  extended  target  or  object  plane  over  which  there  is 
a  continuous  variation  of  contrast.  That  definition  is  directly 
related  to  the  spatial  frequency  bandwidth  associated  with  the 
system  MTF. 
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He  showed  that  the  image  of  a  periodic*  object  plane  pattern 
will  be  spatially  filtered  by  a  system  aperture  of  finite  size, 
so  that  all  Fourier  components  above  a  certain  frequency  in  the 
object  plane  pattern  will  be  removed.  Thus,  there  will  be  a 
minimum  spatial  period  1  such  that  the  image  of  a  sine  wave 
pattern  with  a  period  larger  than  1  may  exhibit  some  modulation, 
but  the  image  of  one  with  a  period  smaller  than  1  will  always 
have  zero  modulation,  i.e.,  constant  intensity  across  the  image 
plane. 

He  showed  that  for  incoherent  illumination 


1 


o 


and  for  coherent  illumination 


(10) 


1 


o 


(11) 


Moreover,  his  demonstration  of  Eqs .  10  and  11  makes  it  clear  that 
the  result  cannot  be  affected  by  any  kind  of  aperture  weighting. 
That  is,  superresolution  techniques  will  not  improve  resolution 
for  extended  objects,  although  apodization  or  other  kinds  of 
processing  may  improve  contrast  or  enhance  selected  spatial 
frequency  components. 

This  result  also  implies  that  the  resolution  limit  for  dis¬ 
tinguishing  details  of  an  extended  target  is  nearly  the  same  as 
that  for  distinguishing  point  targets  without  the  use  of  super- 
resolution  when  the  illumination  is  incoherent  but  is  25% 
poorer  when  the  illumination  is  coherent.  According  to  Eq.  11, 


For  a  finite  field  of  view  this  is  no  restriction  because  the 
pattern  can  always  be  expanded  in  a  Fourier  series  inside  the 
field  of  view. 
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with  a  3.2  mm  wavelength  radar  at  a  distance  of  1  km,  in  order 
to  obtain  the  minimum  spatial  period  of  0.33  m  required  for  tar¬ 
get  identification  it  would  be  necessary  to  have  an  aperture 
at  least  19.^  m  in  diameter. 

f.  Field  Of  View.  The  field  of  view  may  be  defined  as  the 
angle,  measured  from  the  center  of  the  entrance  pupil,  that  is 
subtended  by  the  part  of  the  extended  target  for  which  the  sys¬ 
tem  actually  produces  imagery.  Thus,  the  return  from  any  reso¬ 
lution  element  on  the  target  within  the  field  of  view  must  con¬ 
tain  rays  that  pass  through  the  entrance  pupil  aperture.  A 
target  point  (resolution  element)  may  be  outside  the  field  of 
view  either  because  it  is  not  illuminated  or  because  none  of  the 
rays  in  the  pencil  of  reflected  rays  that  it  generates  pass 
through  the  entrance  pupil  aperture.* 

The  pencil  of  reflected  rays  generated  by  each  target  reso¬ 
lution  element  in  the  field  of  view  produces  a  phase  and  ampli¬ 
tude  distribution  over  the  entrance  pupil  aperture.  Except  for 
aberrations  mentioned  earlier,  the  phase  contributed  to  each 
point  on  the  aperture  is  proportional  to  the  path  length  along 
the  ray  connecting  the  resolution  element  with  the  point. 

The  (coherent)  electromagnetic  field  associated  with  a 
pencil  of  rays  can  be  expressed  as  a  superposition  of  plane 
waves,  each  of  which,  in  accordance  with  its  propagation  direc¬ 
tion,  contributes  a  linear  phase  variation  over  the  aperture. 

The  range  of  spatial  frequencies,  and  therefore  of  propagation 
directions,  of  the  plane  waves  that  contribute  to  the  imagery 
finally  produced  by  the  system,  is  limited  by  the  system  MTF, 
which  ordinarily  vanishes  outside  a  finite  pass  band.  Thus, 
the  MTF  limits  the  system's  field  of  view  as  well  as  its  reso¬ 
lution  and  contrast  modulation. 


Jt 

In  this  connection,  it  should  be  observed  that  the  definition 
of  the  aperture  must  take  into  account  any  of  the  system  stops 
that  may  prevent  rays  from  reaching  the  image  plane. 


22 


g.  Array  Size.  As  shown  in  Appendix  A,  a  system  can  form 
the  image  of  each  target  resolution  element  by  applying  an  appro¬ 
priate  phase  correction  to  the  electromagnetic  field  distributed 
over  the  entrance  pupil  aperture.  It  is  possible  to  carry  out 
this  image  processing  even  if  the  distribution  is  only  given  at 
discrete  sample  points  as  long  as  the  samples  are  dense  enough 
to  satisfy  the  Nyquist  rule:  the  density  must  be  such  that  there 
are  at  least  two  samples  per  cycle  at  the  local  spatial  fre¬ 
quency  corresponding  to  the  aperture  phase  distribution. 

While  for  a  specified  resolution  the  size  of  the  aperture 
covered  by  an  array  of  sensors  is  no  different  from  that  re¬ 
quired  for  continuous  detection,  the  number  of  array  elements 
required  and  their  spacing  depends  upon  the  field  of  view  as 
!  well  as  resolution.  The  total  number  of  elements  needed  and 

their  density  distribution  for'  a  given  resolution  and  field  of 
view  are  calculated  in  Appendix  B.  Because  the  phase  distri¬ 
bution  contributed  by  a  target  resolution  element  over  the 
entrance  pupil  is  such  that  spatial  frequencies  are  lower  at 
the  center  of  the  aperture  than  near  its  perimeter,  the  element 
density  can  be  reduced  at  the  center  without  violating  the 
^  Nyquist  rule.  Thus,  the  total  number  of  elements  needed  is 

smaller  than  would  be  required  to  accomodate  the  maximum  spatial 
frequency  in  the  aperture. 

In  order  to  avoid  multiple  beams,  the  usual  rule  for  spac- 
l  ing  elements  In  an  antenna  array  Is  to  locate  them  no  more  than 

one-half  wavelength  apart.  According  to  Appendix  B,  the  total 
number  of  elements  resulting  from  an  application  of  this  rule 
will  always  be  greater  than  the  number  required  for  Imaging  as 

•  long  as  the  field  of  view  that  must  be  covered  is  less  than  90 
degrees  and  the  distance  to  the  object  plane  is  at  least  as 
large  as  the  diameter  D  of  the  aperture. 

If  Nq  is  the  minimum  number  of  elements  required  in  an 

*  array  for  a  field  of  view  equal  to  e,  and  is  the  number  that 
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would  result  from  elements  spaced  one-half  wavelength  apart, 
then  the  ratio  of  these  numbers  is  given  by 


if  the  wavelength  X  is  3.2  mm.  Thus,  in  this  case 

-k  8 

N0  *  1.43  x  10  x  1.17  x  10  *  16,700  elements 

would  be  sufficient. 
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If  a  uniform  distribution  of  elements  over  the  array  is 

postulated,  then  their  spacing  must  satisfy  the  Nyquist  rule 

at  the  largest  spatial  frequency  in  the  aperture.  In  that  case, 

according  to  Appendix  B,  the  number  N  of  elements  is  related 

m 

to  the  ratio  NQ/N^  by 

2 

Ne  6n2  +  8n  +  3 

\  '  N1  6  (n  ♦  l)5  ' 
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Prom  this  relation  it  follows  that  for  the  example  just  con¬ 
sidered  Nm  ~  26,200  elements,  an  increase  of  57%  over  the  mini¬ 
mum  number  that  is  actually  required. 

h.  SAR  Ambiguity  Limits.  A  synthetic  aperture  radar  may 
produce  range  ambiguities  if  the  pulse  repetition  rate  fp  is  too 
high  and  may  produce  azimuth  ambiguities  if  fp  is  too  low  (cf. 
Ref.  11).  The  azimuth  ambiguities  occur  because  the  sampling 
rate  does  not  meet  the  Nyquist  rule  for  the  highest-resolution- 
element  Doppler  frequencies  that  are  generated  during  the  inte¬ 
gration  time  T. 


The  effect  is  to  limit  the  range  coverage  to  a  maximum 
unambiguous  range  Ru  and  the  platform  velocity  v,  which  deter¬ 
mines  the  time  that  it  takes  to  generate  the  synthetic  antenna 
length  L,  i.e.,  the  integration  time  T.  For  a  system  free  of 
ambiguities,  these  parameters  must  satisfy  the  inequalities 


r  <  c  ,  c  6 
Ru  *  2f“  <  W  ’ 

P 

where  6  is  the  minimum  resolvable  distance  in  azimuth, 
second  inequality  is  equivalent  to 


(14) 

The 


f 


P 


(15) 
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If  the  resolution  limit  6  required  to  identify  a  tank  is 
0.33  m,  then  for  a  platform  velocity  of  100  m/sec,  according  to 
Eq.  15,  f  must  be  at  least  600  Hz,  and  the  unambiguous  range 
Ru  obtained  from  Eq.  14  is  250  m. 

In  principle,  the  largest  possible  platform  velocity  should 
be  used  in  order  to  minimize  the  integration  time,  i.e.,  the 
time  T  that  it  takes  to  generate  the  synthetic  aperture.  Mini¬ 
mizing  T  should  tend  to  minimize  the  phase  error  due  to  target 
motion  and  deviations  of  the  platform  from  its  prescribed  path. 

The  integration  time  for  a  target  resolution  element  at 
range  R  is  given  by 


For  a  target  range  of  1  km  and  a  real  antenna  diameter  of  25  cm 
in  the  example  Just  considered,  the  integration  time  is  there¬ 
fore  0.128  sec.  Therefore,  each  target  resolution  element  will 
be  sampled  77  times.  Since,  for  a  range  of  1  km,  increasing  fp 
by  a  factor  of  250  or  less  will  not  cause  range  ambiguities  to 
develop,  it  is  possible  to  obtain  as  much  as  23  dB  improvement 
in  the  S/N  by  increasing  the  pulse  repetition  rate. 
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IV.  CONCLUSIONS 

A.  CONVENTIONAL  IMAGING  RADAR  CAPABILITY 

Because  of  performance  limitations  that  are  inherent  in 
radar  imaging  systems,  as  a  practical  matter  the  imagery  that 
can  be  obtained  with  a  wavelength  of  3.2  mm  will  not  be  good 
enough  for  identifying  tank-size  targets  at  ranges  even  as  small 
as  1  km,  unless  the  system  used  is  a  synthetic  aperture  radar. 

The  major  impediment  is  the  need  for  an  effective  antenna  aper¬ 
ture  diameter  of  at  least  19.5  m  in  order  to  obtain  the  necessary 
resolution. 

For  the  most  part,  tactical  vehicle  target  surfaces  should 
be  sufficiently  rough  (that  is,  the  number  of  scattering  centers 
should  be  sufficient)  to  make  detailed  imaging  possible  with 
millimeter  waves  at  any  aspect.  For  tanks  and  other  gun-bearing 
vehicles,  however,  gun  barrels  are  apparently  too  smooth  to 
appear  in  the  imagery  at  most  aspects,  except  where  the  viewing 
angle  is  normal  to  the  barrel. 

If  the  radar  antenna  consists  of  a  two-dimensional  array 
of  sensors  over  a  19.5  m  diameter  circular  aperture,  the  minimum 
number  of  array  elements  it  would  require  is  about  16,700. 
However,  this  number  will  be  sufficient  only  if  the  illumination 
can  be  restricted  to  a  field  of  view  covering  not  much  more  than 
a  single  target. 

The  estimate  of  16,700  elements  for  the  19.5  m  array  diam¬ 
eter  is  based  on  the  assumption  that  the  spacing  between  sensors 
Is  nonuniform  and  is  no  smaller  than  is  required  to  accommodate 
the  expected  variation  of  phase  distribution  received  over  the 
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aperture.  For  a  uniform  spacing  sufficient  to  satisfy  the  Nyauist 
rule  for  the  minimum  sampling  density,  the  estimate  would  be 
26,200  elements.  At  least  for  small  fields  of  view,  it  is  always 
possible,  in  principle,  to  reduce  the  number  of  elements  required 
in  an  array  by  using  nonuniform  rather  than  uniform  spacing. 

Superresolution  techniques  exist  for  improving  the  resolu¬ 
tion  of  point  targets  beyond  the  basic  limitation  ordinarily 
imposed  by  the  aperture  size.  However,  such  techniques  fail 
when  applied  to  extended  targets,  for  which  modulation  transfer 
is  a  more  appropriate  measure  of  image  quality  than  resolution 
as  it  is  usually  defined. 

The  difference  between  the  limit  of  resolution  in  the  cases 
of  coherent  and  incoherent  illumination  is  greater  for  extended 
targets  than  it  is  for  point  targets.  Thus,  the  conclusion  of 
V.  Corcoran  in  Ref.  7*  regarding  this  difference  (that  the  reso¬ 
lution  limit  for  the  coherent  case  is  1.41  times  as  large  as  for 
the  incoherent  case)  is  not  valid  for  extended  targets. 

B.  SAR  CAPABILITY  AND  POTENTIAL  PROBLEMS 

In  principle,  synthetic  aperture  radar  can  achieve  the 
image  quality  necessary  for  identifying  military  targets  such 
as  tanks  at  ranges  much  larger  than  1  km.  If  the  targets  are 
stationary,  the  only  essential  limitation  appears  to  be  the  amount 
of  power  available. 

Phase  errors  due  to  unpredicted  target  and  platform  motion 
are  always  a  problem  for  SAR.  Since  error  sensitivity  varies 
inversely  with  wavelength,  a  millimeter-wave  system  would  be  an 
order  of  magnitude  more  sensitive  to  motion  effects  than  would 
an  X-band  system.  However,  for  the  same  resolution  the  milli¬ 
meter-wave  aperture  would  be  smaller,  and  therefore  errors 
would  accumulate  over  a  shorter  time. 


No  classified  material  from  this  classified  reference  has  been 
used  in  this  paper. 
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It  should  be  possible  to  measure  the  target's  gross  trans¬ 
lational  motion  and  make  appropriate  corrections,  but  effects 
due  to  other  types  of  target  motion,  such  as  pitch  and  yaw,  will 
still  be  a  problem.  SAR  systems  ordinarily  provide  some  correc¬ 
tion  for  the  various  errors  due  to  platform  motion. 

The  effect  of  target  and  spurious  platform  motion  on  image 
quality  is  also  a  function  of  the  integration  time,  which  should 
be  significantly  smaller  for  the  contemplated  ranges  (1  to  10  km) 
than  it  is  in  the  usual  X-band  applications.  Increasing  the 
platform  speed  will  help  reduce  the  integration  time  still  further. 
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IMAGE  PROCESSING  FOR  A  STARING  SYSTEM 


A.  PRELIMINARY  CONSIDERATIONS 


The  scattering  of  a  monochromatic  electromagnetic  field 
from  a  target  whose  physical  dimensions  are  large  compared  to 
the  wavelength  results  in  a  spatial  distribution  of  the  form 

E  -  £  A  eik*n  ,  (A-l) 

s  n  n 

in  which  each  subscript  indicates  the  contribution  from  a  single 
scattering  center  on  the  target.  The  coefficients  An  in  (A-l) 
are  complex  vectors  that  vary  slowly  over  spatial  regions,  at 
least  where  the  sensor  that  detects  the  field  is  assumed  to  be 
located.  The  constant  k  is  the  wave  number,  related  to  the 
wavelength  X  by 


and  each  ohase  factor  0  is  a  scalar  function  that  determines 

n 

a  one-parameter  family  of  constant  phase  surfaces,  or  wavefronts, 
given  by 

4>n(x,y,z)  *  constant. 

The  objective  of  an  imaging  system  is  to  separate  out  the 
individual  terms  in  (A-l)  in  order  to  obtain  the  contribution 
of  relative  intensity,  proportional  to  | A  |  ,  and  the  relative 
position  of  each  scattering  center,  which  acts  as  a  source  for 
the  intensity.  The  result,  ideally,  will  disclose  a  uniform 


A-3 
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distribution  of  resolution  elements  that  cover  the  target  image 
and,  by  virtue  of  the  intensity  variation  from  one  to  the  other, 
provide  a  pattern  of  contrast  that  characterizes  the  target  geo¬ 
metrically  . 

The  imaging  process  relies  upon  operations  performed  on  the 

phase  functions  <p  .  For  this  reason  detection  of  the  scattered 

n 

field  must  take  place  over  a  region  that  is  large  enough  to 
accommodate  sufficient  variation  in  the  $  to  make  it  possible 
to  distinguish  one  from  another. 

It  will  be  assumed  that  detection  of  the  scattered  field 
takes  place  within  an  entrance  pupil  in  the  plane  z  =  0.  How¬ 
ever,  the  designation  of  the  entrance  pupil  is  more  or  less 
arbitrary;  It  need  not  be  located  where  the  physical  observation 
of  the  field  actually  occurs  as  long  as  the  amplitude  and  phase 
distribution  over  it  is  a  proper  extrapolation  of  the  observed 
field.  Since  the  sensor  area  must  be  of  finite  extent,  the  field 
extrapolated  onto  the  entrance  pupil  Is  assumed,  in  any  case,  to 
vanish  outside  a  two-dimensional  region  referred  to  as  the 
aperture . 

Huygens'  principle  governs  how  a  field  distribution  over 
the  entrance  pupil  would  determine  the  field  at  a  point  (xq, 
yQ,  zq)  on  some  plane  z  ®  zQ  within  the  Imaging  system.  The 
contribution  of  a  single  term  in  (A-l)  would  produce  a  term  of 
the  form 


EQ  =  fJ'Ae^^^'dxdy  ,  (A-2) 

aperture 


V(x- 


+  <y-yo>2  +  zo2 


A-4 


where 


For  most  positions  of  the  ooint  (x  ,  y  ,  z  )  the  anolitude 
effect  of  the  variation  of  the  phase  k(l+Ft)  in  (A-2)  is  assumed 
to  be  large  compared  with  the  variation  in  A  over  the  aperture. 
This  phase  variation,  in  fact,  causes  the  integrand  to  oscillate 
many  times  over  the  integration  region  so  that  its  contributions 
to  the  value  of  the  integral  tend  to  cancel  each  other.  For 
such  positions  of  (xq,  yQ,  zq)  the  intensity  |E  |  given  by  (A-2) 
will,  therefore,  be  relatively  small. 

However,  if  for  other  positions  of  (x  ,  y  ,  z  ),  e.g.,  in 

A  A  A  ^  ^  ^ 

a  small  region  about  some  point  (x,  y,  z),  the  phase  k(^+R)  is 
constant  or  nearly  constant  over  the  aperture,  all  contributions 
of  the  integrand  to  the  value  of  the  integral  tend  to  add  in 
phase.  In  this  case  the  intensity  will  be  relatively  large,  and 

AAA 

the  neighborhood  of  the  point  (x,  y,  z)  may  be  regarded  as  a 
focal  region. 

The  object  of  image  processing  is  to  modify  the  phase  over 
the  entrance  pupil  In  such  a  way  that  a  separate  and  distinct 
focal  region  is  created  for  each  term  In  (A-l).  That  Is,  each 
4>n  Is  modified  for  this  purpose  by  the  addition  of  a  phase  shift, 
a  process  which  Is  accomplished  in  classical  optical  systems  by 
means  of  a  lens. 

B.  PLANE  WAVE  IMAGING  (DISCRIMINATING  ANGLES  OF  ARRIVAL) 

The  wavefront  <t>  *  constant,  radiated  from  a  single  scatter¬ 
ing  center,  will  be  approximately  a  plane  wave,  i.e.,  will  have 
the  form 

4>(x,  y)  *  px  +  qy  (A— 3 ) 

over  the  aperture  if  the  distance  to  the  scattering  center  is 
sufficiently  large  compared  to  the  aperture  diameter  D.*  In 
( A— 3 )  p  and  q  are  direction  cosines  of  the  propagation  direction 

•  2D^ 

The  distance  should  be  of  the  order  -r—  or  more,  which  guarantees 

that  the  error  in  *  across  the  aperture  Is  no  greater  than  1  . 


I 

k 


and,  as  such,  provide  the  angle  of  arrival  of  the  wave.  That 
is,  the  wave  propagates  in  the  direction  of  the  unit  vector  with 
components  (p,  a,  r),  where 


r  =  V  l-p2-q^  .  (A~4) 

The  definition  (A-4)  implies  that 

P2  +  q2  *  1  .  (A-5) 

Plane  wave  images,  after  processing,  will  appear  in  the 
system’s  focal  plane  located  at 


where  f  is  the  focal  length.  It  will  be  assumed  that  the  focal 
plane  is  in  the  Fresnel  diffraction  region  of  the  entrance  pupil, 
i.e.,  f  is  large  enough  compared  to  the  diameter  D  of  the  aper¬ 
ture  so  that  the  approximation 


R  ~ 


2  2 

x  +  y 

2f 


(XQ  *  *  yQ  r) 


x  2  +  y  2 
o  Jo 

2f 


+  f  (A-6 ) 


is  valid  In  (A-2).  This  Is  equivalent  to  requiring  that  the 
F-stop,  which  Is  the  ratio  of  the  focal  length  to  the  aperture 
diameter,  be  large.* 


The  image  processing  for  plane  waves  consists  of  modifying 
the  phase  from  <t> ,  given  by  (A-3),  to 

♦  •  ♦  -  X-2f -y~  =  px  +  qy  -  —  ~  2  ■  .  (A— 7 ) 

Then  the  total  phase  in  the  Integrand  of  (A-2)  becomes 

*  +  R  -  (p  -  +  (q  -  ^)y  +  -Tf-y-  +  f  (A-8) 


because  of  (A-6). 


« 

According  to  the  classical  Rayleigh  criterion,  the  maximum 
error  in  (A-6)  across  the  aperture  should  be  less  than  K  . 


It  follows  from  (A-8)  that  if  x  ,  yQ  and  zq  are  coordinates 


of  a  point  in  the  focal  plane,  defined  by 

xo  =  fp,  yQ  -  fq,  z0  =  f  » 


(A-9) 


the  phase  ^+R  is  constant  over  the  aperture.  Thus,  for  each  p 
and  q  that  define  an  angle  of  arrival  there  is  an  image  point, 
given  by  (A-9),  where  the  intensity  appears  to  be  focused. 

C.  POINT  IMAGING 

When  a  target  scattering  center  is  a  point  located  at  (x^ 

y  ,  z^)  and  is  not  far  enough  from  the  entrance  pupil  to  be  in 

the  Fraunhofer  region  (for  which  the  distance  would  be  of  the 
2 

order  2D  or  more),  then 
X 


(x-x^)2  +  (y-j^;  + 


( A-10 ) 


If  the  distance  is  far  enough  for  the  target  to  be  in  the  Fresnel 
region,  however,  then  an  approximation  similar  to  (A-6)  is  valid, 
so  that 


x2  ♦ 


X;,x  +  y,y  x2  +  y2 

— - —  +  - - -  +  z,  . 

z^  2z^  1 


(A-ll) 


In  this  case  the  image  plane  will  not  be  the  focal  plane 
and,  therefore,  f  in  (A-6)  should  be  replaced  by  the  more  general 
z  .  Nevertheless,  the  imaging  process  consists  of  applying  the 
same  phase  correction  used  in  (A-7)  for  the  plane  waves;  i.e.. 


+  R  *  $  - 


x2  + 


(A-12) 


Because  of  (A-ll)  and  the  new  approximation  for  R  the  total 
corrected  phase  is  approximately  given  by 


where 


(x2  +  v^) 

+  R  -  a  — — 2  *  +  six  +  s2y  +  y 


(A-13) 


a  ■  —  +  — 

Z1  zo 
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(A-1H) 


(A-15) 


and 


(A-16) 


The  coefficient  a  given  by  (A-14)  will  vanish  if  the  image  plane 
is  identified  with  the  conjugate  plane  determined  by  the  lens 
law  of  classical  optics;  i.e.,  the  value  of  zq  will  be  determined 
by 


1 


Z1 


+ 


0  . 


(A-17 ) 


The  phase  distribution  over  the  entrance  pupil  will  then 
have  the  same  form,  i.e.,  linear,  as  in  the  plane  wave  case,  for 
which  the  distribution  is  given  by  (A-8).  In  the  present  case, 
however,  the  coordinates  of  the  image  point  are  determined  by 
setting  and  82,  which  are  defined  by  (A-15),  equal  to  zero. 
The  result  is 

xo  =  Mxl  »  yo  =  Myl  »  (A-18) 

where,  because  of  (A-17), 

M  =  ,  -  ■■  ■  .  ( A-19 ) 

i  -  z1 

The  quantity  M  is  obviously  the  magnification  factor  for  the 
image . 

D.  RESOLUTION 


For  both  the  plane  wave  and  the  point  image  cases  the  field 
in  the  appropriate  image  plane  is  given  by  a  form  of  (A-2)  in 
which  the  total  phase  of  the  integrand  is  a  linear  function  over 
the  aperture  because  of  (A-8)  on  the  one  hand  and  (A-13)  on  the 
other.  If  the  slowly  varying  amplitude  factor  in  (A-2)  is 


A-8 


assumed  to  be  constant,  then  the  field  in  either  case  will  be 
proportional  to  a  quantity  I  of  the  form 

i  =  f yelk(5x-ny>dx(Jy  _  (A.20 

aperture 

If  the  aperture  is  a  circular  region  of  radius  a,  then 


I 


a  air 


// 


eikPrcos(e-e0)rd6dr 


( A-21 ) 


where 


.  •  . 

The  evaluation  of  the  integral  in  (A-21)  is  well  known  (Ref.  A-l) 

3  Ocap) 

1  -  27ra  •  (A'22) 

where  is  the  Bessel  function  of  order  one. 

The  Rayleigh  criterion  for  the  limit  of  resolution  of  two 
point  targets  is  defined  as  the  separation  for  which  the  peak 
of  the  image  plane  pattern  given  by  (A-22)  for  one  is  at  the 
first  null  of  the  image  plane  pattern  for  the  other.  Since  the 
peak  value  of  I  occurs  where  p  is  zero,  the  smallest  solution 
pQ  of  the  equation 

J1(kapQ)  -  0  (A-23) 

Is  the  limit  of  resolution  according  to  the  Rayleigh  criterion. 
That  solution  Is  the  well-known  result 

p  .  »  X-‘2.2X-  .  (A-24) 


A-9 


E.  DEPTH  OF  FIELD 


The  image  of  a  point  target  is  focused  in  the  conjugate 
plane  determined  by  the  lens  law  of  (A-17).  That  is,  the  total 
phase  i^+R  of  the  integrand  in  (A-2)  is  constant  over  the  aper¬ 
ture  as  long  as  xq  and  yQ  satisfy  (A-18)  and  zq  satisfies  the 
lens  law  (A-17).  The  image  will  effectively  remain  in  focus 
even  if  zq  deviates  from  its  conjugate  plane  value  by  an  amount 
Az  as  long  as  the  corresponding  deviation  in  the  value  of  the 
phase  factor  \|<+R  remains  less  than  the  Rayleigh  error  limit  of 
r-  .  The  interval  determined  by  the  excursions  Az  that  do  not 
violate  the  Rayleigh  limit  can  be  defined  as  the  depth  of  focus. 
If  the  image  plane  at  z  =  zQ  remains  fixed  and  Az  is  the  excur¬ 
sion  of  the  target  from  its  longitudinal  position  at  z^  where 
the  object  plane  is  located,  the  corresponding  interval  for  which 
the  Rayleigh  limit  is  not  exceeded  can  be  defined  as  the  depth 
of  field. 

Accordingly,  if  B1  and  B2  are  set  equal  to  zero  and  z^  Is 
replaced  by  z1+Az  in  (A-17),  then,  except  for  a  term  that  is 
constant  over  the  aperture. 


<|>  +  R  -  £  (x2  +  y2)  , 


(A-25) 


where 


Az 

z1(z1  +  Az)  * 


(A-26) 


The  maximum  phase  change  occurs  at  the  aperture  edge  where 


and,  therefore,  is  given  by 


ip  +  R 
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At  the  Ravleieh  Unit  this  would  be  such  that 


1.6., 


U  a2  .  +  X  . 
7  a  "  4  H  » 


±  _±_  *  4  2X  < 
2a2  D2 


(A-27) 


The  depth  of  field  interval  is,  therefore,  defined  by 


2  2 

2XzJ  2Az( 

Z1  ~  ~2  s  R  i  +  — 2 

x  D  +  2Az1  1  D  -  2Xzx 


(A-28) 


where  R  is  the  distance  to  the  target  plane  from  the  entrance 
pupil  plane  at  z  *  0  after  the  target  has  moved.  If 


D  >>  2Xz,  , 


then  this  becomes 


2Az2  2Az? 

Z1  “  2~  s  R  s  zi  +  ~ 

D 


(A-29) 


The  image  point  field  given  by  (A-2)  with  (A-25)  is  pro¬ 
portional  to  the  integral  I  given  by 


2tt  a 


l  l' 


-ik£r 


rdrd8  , 


(A-30) 


if  the  aperture  is  circular  and  the  slowly  varying  amplitude  is 
assumed  to  be  constant.  The  integral  In  (A-30)  is  easily 
evaluated: 

iffW*2  2 

t _ 2X  _  tA  - j  us  /»  n  \ 


1  ~  T 


sm^f-. 


(A-31) 


A-ll 


At  its  peak,  where  y  =  0 , 


the  intensity  is  proportional  to 


1 1 1 2  =  ttV  .  (A-32) 

At  the  extreme  edge  of  the  depth  of  field  interval,  where  u  has 
the  value  given  by  (A-27),  the  intensity  is  proportional  to 

1 1 1 2  =  8a4  .  (A-33 ) 

The  ratio  of  the  intensity  corresponding  to  (A-33)  to  the  peak 
intensity,  corresponding  to  (A-32),  is  8/n  or  0.81.  This 
corresponds  to  a  drop  of  19%  for  the  intensity  from  its  value 
where  the  target  is  at  the  position  of  best  focus  to  its  value 
when  the  target  is  at  the  edge  of  the  depth  of  field  interval. 

F.  REALISTIC  SCATTERING  CENTER  IMAGERY 

The  image  processing  discussed  in  this  appendix  thus  far  is 
based  on  two  assumptions:  (1)  the  scattering  center  is  a  point; 
(2)  the  F-stop  is  large  enough  for  the  rays  that  enter  the  system 
to  be  regarded  as  paraxial.  In  practice  both  assumptions  will 
be  incorrect  to  some  degree,  and  to  the  extent  that  they  are 
incorrect  both  lead  to  phase  errors,  or  optical  aberrations, 
that  cause  the  image  plane  diffraction  patterns  to  be  more  com¬ 
plicated  than  predicted.  The  effect  is  to  degrade  the  resolu¬ 
tion  and  make  the  images  appear  to  be  out  of  focus. 

The  aberrations  due  to  violation  of  the  second  assumption 
are  predictable  and  well  known  in  classical  optics.  As  past 
experience  has  demonstrated,  it  is  possible  to  correct  them  when 
necessary. * 


1 - 

A  general  approach  that  has  often  been  useful  in  synthetic 
aperture  radar  processing  is  to  regard  the  correction  of 
predictable  aberrations  as  a  problem  of  designing  a  spatial 
two-dimensional  (or  three-dimensional)  matched  filter. 


A-12 


The  aberrations  due  to  the  first  assumption  are  most  trouble¬ 
some  when  the  illumination  is  coherent.  They  are  apparently  the 
cause  of  the  well-known  speckle  effect,  which  is  averaged  out 
by  incoherent  illumination.  Since  they  are  not  predictable, 
there  appears  to  be  no  way  of  eliminating  them  except  by  some 
kind  of  averaging  process,  such  as  dithering  the  illumination 
source  or  possibly  the  image  plane  position  and  orientation. 
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APPENDIX  3 


STARING  SYSTEM  ARRAY  THINNING 


If  the  spacing  d  between  the  elements  of  an  antenna  array 
is  larger  than  half  the  wavelength  X,  then  the  detection  of  the 
angle  of  arrival  of  a  plane  wave  will  be  ambiguous  unless  it  is 
known  that  the  field  of  view  is  suitably  restricted  to  an  angle 
less  than  180  degrees.  If  it  is  known  that  the  field  of  view  is 
no  larger  than  an  angle  9,  then  the  detection  will  be  unambiguous 
as  long  as 

d  s  — * — =.  .  (B-l ) 

2  sin^ 


From  another  point  of  view,  (B-l)  is  equivalent  to  the 
Nyquist  condition  for  the  minimum  sampling  density  required  in 
the  array  plane  in  order  to  reconstruct  the  incident  wave.  That 
is,  limiting  the  field  of  view  to  the  angle  0  is  equivalent  to 
limiting  the  spatial  frequency  pass  band  of  the  signal  distri¬ 
bution  over  the  array  plane  to  frequencies  bounded  by 


f 


max  55  X  sln2 


(B-2) 


Therefore  (B-l)  is  Just  the  requirement  that  the  sampling  den¬ 
sity  be  at  least  two  samples  per  cycle  at  the  maximum  spatial 
frequency  to  be  detected  by  the  array  elements. 

More  generally,  for  any  coherent  incident  wave,  if  the 
phase  distribution  over  the  array  plane  is  a  function  4>(x,  y), 
then  the  spatial  frequency  at  a  point  (x,  y)  in  the  direction 


B-3 


of  a  unit  vector  £  with  components  (cos<{>,  sinb)  is  given  in 
terns  of  the  two-dimensional  gradient  7$  by 

■  v*  •  (B-3> 

For  an  incident  wave  of  the  form 

E  =  A  elk^(x*  y’  z) 

the  phase  distribution  over  the  array  plane  is  given  by 

$(x,  y)  =  ’Kx>  y>  °)  »  (B-4) 

and  therefore,  according  to  (B-3)» 

fs(x,  y)  =  ^[i|ix(x,  y,  0  )cos<j>  +  i|»y(x,  y,  0)sin<()]  .  (B-5) 

For  example,  in  the  case  of  a  plane  wave 

y,  z)  =  px  +  qy  +  rz  ,  (B-6) 

where  p,  a,  and  r  are  the  direction  cosines  of  the  wave’s 
propagation  direction.  Then,  according  to  (B-5), 

fs(x,  y)  *  ^  (pcos<|)  +  qsin<j>)  *  j-  sin0Q  ,  (B-7) 

where  0  is  the  angle  of  incidence.  Since 

o 


if  the  incidence  direction  is  the  maximum  permitted  by  the 
field  of  view,  (B-7)  is  consistent  with  (B-2). 

If  a  target  point  located  at  (x^,  y^,  z^)  acts  as  a  source 
for  the  incident  wave,  then 

(x ,  y,  2)  -  y  (x-xx)^  +  (y-y1)2  +  (z-zx)2 


B-4 


( B—  8 ) 
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If  r  is  chosen  to  have  the  direction  of  V$,  then  f  will 
»o  s 

have  its  maximum  value  and,  according  to  (B-3)  and  (E-4),  will 
be  Kiven  by 

fs  =  J|7^(x,  y,  0) |  .  ( B— 9 ) 

From  (B-8)  and  (B-9)  it  follows  that 


1  (x-x]_)2  +  (y-y1)2 

*2  [(x-x1)2  +  (y-y.^)2  +  z12] 


which  can  also  be  written 


(B-10 ) 


2  _  1 
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(x-x1)2  +  (y-y1)2 


+  Z1  . 


(B-ll ) 


Suppose  that  the  array  aperture  is  circular.  From  (B-ll) 

p 

it  is  evident  that  the  maximum  value  of  f  on  any  circle  of 

s 

radius  p  about  the  center  of  the  aperture  is  given  by 


where 


f2(p) 

max 


( p+p-j^ ) ' 


[(p+P1)2  +  Zl2] 


(B-12 ) 


P1 


This  is  shown  in  Fig.  B-l. 

Now  the  Nyquist  rule  that  there  must  be  at  least  two  sam¬ 
ples  per  cycle  of  spatial  frequency  in  any  direction  will  be 
satisfied  if  the  area  density  o  of  elements  in  the  array  satis¬ 
fies  the  relation 


o  »  Df 
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max 


(B-13) 
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on  every  infinitesimally  thin  annulus  around  the  center  of  the 

aperture.  On  the  other  hand,  the  total  number  N  of  elements 

5  e 

in  the  array  may  be  defined  as  the  smallest  integer  such  that 


//■ 


dxdy 


aperture 


(B-U) 


Then  it  follows  from  (B-13)  and  (B-l^)  that  the  Nyquist  rule  can 
be  satisfied  if 


Ne  =  4  dxdy  +  n  ’ 


(B-15) 


where 


aperture 


0  s  n  <  1  . 


For  a  circular  aperture  of  radius  a,  after  substituting 
from  (B-12 )  into  (B-15),  it  is  found  that 
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/pdp  -  z2  f - ^ 

J  1 J  ( P+PT)2  +  zf 

L  n  n  ±  -L 


+  h  . 


That  is,  after  evaluating  the  integrals, 

B-7 
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(B-16) 


Ne  =  ^  [a2  +  z22  U(a ,  +  2plZl  V(a,  z^]  +  n  , 


where 


U(a ,  z  ^ )  *  log  I  1  + 


-  log 


(a  +  p1)‘ 


i  a  +  p.  ,  p, 

V(a,  z,  )  =  tan  - ±-  -  tan"-1  -i 

l  zx  Zi 

Prom  the  fact  that  for  small  e 

2  3 

log  (1  +  e) ~  e  -  |—  +  |- 


-1  £ 
tan  e  ~  e  —  +  ... 


it  follows  from  (B-16)  that  for  large  z^,  i.e.. 


z1  >>  a  +  p1  , 


N  -  ,.?q3 

6  (xZl)2 


/ 6n2  +  8n  +  3\ 

\  3  /  ’ 


(B-17 ) 


where 


For  elements  spaced  one-half  of  a  wavelength  apart  the 

2 

linear  density  would  be  j  and,  therefore,  the  area  density 


of  elements 


2 

would  be  4/A  .  In  that  case  the  total  number  N^/2 
over  the  aperture  would  satisfy 


Therefore,  combining  (B-17)  and  (B-l8), 


The  ratio  is  less  than  one  as  long  as 


(B-18) 


(B-19) 


(B-20 ) 


which  is  equivalent  to  the  field  of  view  being  less  than  90  degrees. 

Another  way  to  satisfy  the  Nyquist  rule  would  be  to  use  a 
uniform  density  of  elements  over  the  aperture  with  the  density 
equal  to  what  would  be  required  for  the  largest  spatial  frequency 


This  estimate  is  conservatively  small  because  the  density 
should  be  somewhat  higher  in  order  to  account  for  the  wider 
spacing  that  would  occur  in  the  diagonal  direction. 


B-9 


observed  by  any  part  of  the  array.  According  to  (B-12),  the 
largest  spatial  frequency  is  given  by 


a  +  p, 

f£il  =  ~F=  =  == 

W(a  +  p,)2  +  zf 

(B -22) 

_  (n  +  l)a _ 

xV(n  +  l)2a2  +  z2 


after  a  substitution  based  on  (B-20).  The  corresponding  element 

density  is  given  by  4f2(a),  and  therefore  the  total  number  N 

max  m 

of  elements  in  the  array  is  given  by 


because  of 
that 


(B-22 ) 


Prom  (B-17) 


4ffa4(n  +  l)2 
[  (n  +  l)2a2  +  z2] 

,  (B-19),  and  (B-23), 


(B-23) 

it  follows 


Ne 

% 


6n2  +  8n  +  3 
6(n  +  l)2 


(B-24) 


B-10 


END 

DATE 

FILMED 

1-82 
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